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Abstract. The generalized Jacobi equation is a differential equation 
in local coordinates that describes the behavior of infinitesimally close 
geodesies with an arbitrary relative velocity. In this note we study some 
transformation properties for solutions to this equation. We prove two 
results. First, under any affine coordinate changes we show that the ten- 
sor transformation rule maps solutions to solutions. As a consequence, 
the generalized Jacobi equation is a tensor equation when restricted to 
suitable Fermi coordinate systems along a geodesic. Second, in dimen- 
sions n > 3, we explicitly show that the transformation rule does not in 
general preserve solutions to the generalized Jacobi equation. 



1. Introduction 

Suppose M is a manifold with an affine and torsion free connection V. In 
this setting, the Jacobi equation is a fundamental equation that describes the 
qualitative behavior of infinitesimally close geodesies on M. See for example, 
[LC27j and [Hic65] . One way to derive the Jacobi equation is to consider the 
deviation = x fi (s) — X^(s) between two neighboring geodesic X and x. 
Then the Jacobi equation follows by subtracting the geodesic equations for 
X and x and assuming that P(s) and £*(s) = ^(s) are both infinitesimal 
quantities. This is a standard argument. See for example [Per08| and refer- 
ences therein. The generalized Jacobi equation is derived in the same way, 
but under the weaker hypothesis that only is an infinitesimal quantity. 
That is, the geodesies are assumed to be infinitesimally close, but their rel- 
ative velocity £^ does not need to be small. Under this weaker hypothesis, 
equations for the displacement £ M was derived first in the Lorentzian case 
by Hodgkinson [Hod72] and independently by Ciufolini [Ciu86]. This gener- 
alization of the Jacobi equation have been investigated by several authors, 
specially with applications on astrophysics and cosmology. See for instance 
[Mas75llCM02llPer08j . 

We will work in the setting of an affine and torsion-free connection. In this 
setting, the generalized Jacobi equation was derived by Perlick |Per08| . In 
more detail, a collection of n functions ^ : I — > R along a geodesic X : / — > 
M is a solution to the generalized Jacobi field in coordinates {x m }™Zq if 

(i) ^+rz p (2ex»+ec) + ^e(x p +£ p ) (>+r) = o, 
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where dots indicate ordinary derivatives with respect to the parameter s of 
the central geodesic X^(s) and Tu P are the connection coefficients of V. 

Let us first observe that in contrast to the usual Jacobi equation, the gener- 
alized Jacobi equation is a non-linear equation in the unknown functions 
This makes the analysis of the solution space more difficult. For example, 
the solution space need not be a vector space and in general, there are no 
(known) results that associate a solution £^ to a geodesic variation as for 
the Jacobi equation. 

The generalized Jacobi equation (JTJ) is an equation in local coordinates. If 
V is flat, there exits coordinates where T„ p = [Shi07t Proposition 1.1]. 
Then the generalized Jacobi equation and the usual Jacobi equation both 
simplify into = 0. More generally, Perlick |Per08j has proven that if X 
is a lightlike geodesic in a special class of Lorentz metrics, or planewave 
metrics, then there are coordinates around X where functions £ M satisfy 
equation (pQ) if and only if functions £ M satisfy the usual Jacobi equation. 
When this is the case, it implies that by a suitable choice of coordinates, the 
nonlinear generalized Jacobi equation can be replaced by the linear Jacobi 
equation. Since the Jacobi equation transforms tensorially, this motivates 
a further understanding of the transformation properties of the generalized 
Jacobi equation. Namely: 

Suppose X : I — > M is a geodesic in overlapping coordinates x^ and 
x^. If functions I — > R solve equation (pQ) in coordinates x^, 
does there exist a transformation rule ^ h-> £ m such that functions 
£ M solve equation ([1]) in coordinates x^l 

In this note we will study this question for the tensorial transformation 
rule £^ = ffi7^- We will prove two results. First, in proposition 14.11 we 
show that the tensorial transformation rule for £ l preserves solutions to 
equation (QJj for affine coordinate changes. As a consequence, equation ([I]) 
is a tensorial equation when restricted to suitable Fermi coordinate systems 
along a geodesic (see Proposition 12. 2[) , This motivates the use of Fermi 
coordinates for the study of the generalized Jacobi equation as in [ICM02, 
IPer08| . Second, in proposition 14.21 we explicitly show that in dimensions 
n > 3, the tensorial transformation rule for £^ does not in general preserve 
solutions to equation ([1]). Thus, if there exists a transformation rule £^ — > £^ 
for solutions to equation (pQ), it is not the tensorial transformation rule. 

The organization of this note is as follows. In Section [2] we review the 
necessary theory for affine connections and Fermi coordinates. In Section 
[3] we summarize the derivation of the Jacobi equation and the generalized 
Jacobi equation |LC27j . Lastly, in Section U] we prove proposition 14.11 and 
proposition 14.21 described above. 

2. Preliminaries 

Let M be a smooth manifold of dimension n > 2. By TM we denote the 
tangent bundle with projection tt: TM — > M. The tangent space at a point 
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p E M is defined as T p M = ir 1 (p). Throughout the paper we denote by I 
an interval in R. We will also use the Einstein summing convention. 

We assume that M is endowed with an affine connection V. Thus, in each 
coordinate chart (U, V is represented by connection coefficients Yy a and 
if Ty a and Ty a represent V on overlapping coordinates x^ and x^, we have 
transformation rules 



dx x ~x dx a dx@ d 2 



x 



v ; dx a va a(3 dx v dx° dx v dx° ' 

A connection is torsion- free if Ty a = T% u . If A is a curve X : I — > M, then 
X is an (affinely parameterized) geodesic if (i) X is a regular curve, that is, 
the tangent X is never zero, and (mJ in each coordinate chart (£7, x M ) that 
overlaps X we have 

(3) A^ + r^A^jr = 0, 
where (A^(s)) are components for A in coordinates x^ . 

2.1. Fermi coordinates. In this section we collect some results on Fermi 
coordinates. These are local coordinates for a tubular neighborhood around 
a geodesic [MM63j . Suppose X C R is an open interval and A: X — >• M 
is a geodesic of an affine connection on M, and suppose II is an (n — 1)- 
dimensional vector space in TW So )M for some so £ I such that II is com- 
plementary to A(so), that is, 

(4) T x(so) M = span{A(s )} II. 

Let {ei, . . . ,e n _i} be a basis for EL By parallel transport we can extend 
each vector ei £ Tjf( S0 )M into a vector field ej(s) along A. Since parallel 
transport is a linear isomorphism, it follows that 

(5) T x{s) M = span{A(s),ei(s), . . . ,e„_i(s)} 
for all s £ I. For s £ I, let / be the map 

fn-l 

(6) /(a,* 1 ,...,*"" 1 ) = exp \J2 



- r; 

. 1=1 



defined for z , . . . , z n € R for which the right hand side is defined. 

Proposition 2.1. Suppose X is an open interval, A: I — > M is a geo- 
desic for an affine connection V, and {e^Zi an d f are as * n equation (0). 
Moreover, suppose X has no self-intersections, and I is a proper open subset 
I C X such that I is compact. Then there exists an open neighborhood of the 
origin B C R n_1 such that f restricts to a diffeomorphism f: I X B M 
onto its range. 

Proof. Since I is compact, we can find an open ball B C R n_1 containing 
such that / : I x B — > M is smooth. For s E / we have 

!^( S ,0,..., 0) = A( S ), 
df 

^-(s,0,...,0) = ei (s), i€{l,...,n-l}. 
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Since {X(s), ei(s), . . . , e n _i(s)} are linearly independent for all s £ I, the 
inverse function theorem implies (after possibly shrinking B) that f: I X 
B — > M is a local diffeomorphism onto its range. The result follows by 
|5pl79] p. 345, Lemma 19]. □ 



When (s, z , . . . , z n ~ l ) £ IxB are as in Proposition ^. ll we say that (s, z , . . . , z n ~ l ) 
are Fermi coordinates along X : I — > M. These coordinates are determined 
by the geodesic X : I — > M, the initial point X(sq) and the set of vectors 
{e^so)}™!] 1 in equation ([5]). One can prove that when the connection is 
torsion-free, then all Christoffel symbols T^p vanish on the central geodesic 
in Fermi coordinates. The next proposition shows that if two Fermi coordi- 
nates systems are determined by the same initial complementary hyperplane 
IT, then the Fermi coordinates differ by an affine coordinate transformation. 
Here, two overlapping coordinates (U, x^) and {U ,x^) are related by an 
affine coordinate transformation if 

(7) x^(x ,...,x n ^) = A^ + C 
for some constants (A^)"~L and (C^)™^. 

Proposition 2.2. Suppose V is an affine connection on M and X : I — )■ M 

and X : / —> M are two geodesies that differ by a reparameterization. More- 
over, suppose (s, z 1 , . . . , z n_1 ) and (s, z 1 , . . . , i™ -1 ) are Fermi coordinates 
along X and X that correspond to the same initial complementary hyper- 
plane II. Then Fermi coordinates (s,z , . . . , z n ) and (s,z , . . . , i™ -1 ) are 
related by an affine coordinate transformation on their common domain. 
More precisely, there are constants A, B and an invertible matrix (T^)™~ =1 
such that 

n-1 

(8) s = As + B, F = ^T^/or/iG{l,...,n-l}. 

u=l 

Proof. Let / and / be maps /: I X B — > M and / : / x B — > M as in 
equation Q that define the two Fermi coordinates. Let ip: I — > I be the 
reparameterization, so that Xotp = X. Writing out the geodesic equation for 
X and Xotp shows that tp(s) = As + B for some A £ R\{0} and B £ R. By 
assumption, there are so £ I and So G I such that vectors {j^i\x(s )}p]Zi 
and {g^rlx(j )}^=i span the same hyperplane. Thus X(sq) = X(So) and 
there exist an invertible matrix T = (T^)™"^ G R( n - 1 ) x ( n - 1 ) such that 



< 9 > h 



n-1 



X(s) „ =1 ° Z 



, n £ {1,... ,n- 1} 

X(s) 

at s = so. Since both sides in this equation are parallel vectors along X, 
and since equality holds for s = so, it follows that equation ([9]) holds for all 
s£l. 

Suppose (s, z , . . . , z n_1 ) and (s, z , . . . , z n ~ l ) are Fermi coordinates for the 
same point, so that 

(10) f{s,z\...,z n - 1 ) = /(a, i 1 ,...,?- 1 ). 
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Then equations © and ([9|) imply that 

(n—l n— 1 

As + B^Ty,...,^^- 1 ; 
jU=l jU=l 

Since / is a bijection onto its range, equations (fT0|) - (fTT|) imply that equation 
© holds. □ 

Let us consider the case when (M, g) is a pseudo-Riemann manifold of index 
1 and V is the Levi-Civita connection of g. If X(sq) is timelike or spacelike, 
then a suitable hyperplane IT is given by IT = (X(so)) 1 - |Q'N831 p. 49]. In 
this case, the metric can further be made diagonal along the geodesic in 
Fermi coordinates |LC27j . On the other hand, if X(sq) is lightlike, we have 
X(sq) £ (X(so)) ± and the choice II = (X(sq)) 1 - is not possible. 

3. Equations of geodesic deviation 

In this section we describe three equations for the behavior of nearby geodesies: 
the exact deviation equation, the Jacobi equation, and the generalized Jacobi 
equation. Throughout this section we assume that V is an affine and torsion- 
free connection on M. 



3.1. The exact geodesic deviation equation. Suppose x: I —> M and 
X : I — > M are two geodesies that are contained in one coordinate chart 
(U,x»). If locally x(s) = (x»(s)) and X(s) = (X>"(s)) for s £ I, let £:/—)• 
W 1 be the displacement between the geodesies defined as 

(12) e(s) = x^(s)-X^{s), s€l. 

Since x and X are solutions to the geodesic equation, it follows that 

(13) e + r^(x + o (x v + r) (x° + f 7 ) - r^(x) X a X v = 0. 

We will refer to (fl~3j) as the exact geodesic deviation equation (see [Per08j). 
Equation ()13p is an exact equation in the sense that its derivation does not 
involve any approximations. However, the geometric analysis of equations 
(fl~2|) - (fT3j) becomes difficult since both equations are defined in local coordi- 
nates and, moreover, both equations involve two points on the manifold. 

That components £^ defined by equation (|12p do not define a vector field 
along X : / — > M can be seen as follows. If 1R 2 is equipped with the Euclidean 
metric and Cartesian coordinates (x l ,x 2 ), then curves X(s) = (s,0) and 
x(s) = (s, 1) for s > are geodesies and functions £ M : / — > M 2 are given by 
£(s) = (0, 1). However, in polar coordinates x 1 = r, x 2 = 9, the definition of 

£ yields eO) = (x/TT^-^tan- 1 ^)). It follows that ? (s) = (g) \ x{s) & 

is not satisfied for all s > 0, and functions £ M in equation (|12p do not in 
general transformation as a tensor. 

Even if functions £**(s) in equation (j 12 j) do not transform as a tensor in 
general, it turns out that if we restrict to suitable Fermi coordinates along 
a geodesic, then functions £°, . . . ,£ n_1 transform as a vector. To see this, 
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suppose X : I — > M is a geodesic and (s, z , . . . , z n_1 ) are Fermi coordinates 
along X as in proposition 12.11 If x : I — >• M is another geodesic that can 
be written as x(s) = (s, z 1 (s), . . . , z n ~ l (s)) in these Fermi coordinates, then 
functions £°, . . . , in equation (fT2|) are given by 



(14) e(s) 




for = 0, 

for fj, € {1, . . . , re — 1}. 

If (s, i 1 , . . . , S 11-1 ) are other Fermi coordinates along X as in proposition 
and £^ are defined by equation (fT2|) in these coordinates, then equations ([6]) 
and ([8]) show that ^ and £ M transform as a vector along X. 

Since we have not made any approximations when deriving equation (|13j) . 
it turns out that we can transform solutions from one coordinate system to 
another. However, this will lead to a highly non-standard transformation 
rule. To see this, suppose X ll (s) and x^(s) are geodesies as above in coordi- 
nates x M whence £ M (s) in equation (fT2"j) solve equation (fTBj) . Thus, if x M are 
overlapping coordinates and X^(s) and x^{s) represent the same geodesies 
in these coordinates, then functions £^(s) = x^{s) — X^{s) solve equation 
(|13|) in ^-coordinates. Then 

(15) e{s) = {xox^f (Z\s)+X x (s))-X»(s), 

and the above equation gives a transformation rule £ M i— > £^ for solutions to 
the exact deviation equation along a geodesic X : I — > M. 

3.2. The Jacobi equation. The usual approach to analyze the qualitative 
behavior of nearby geodesies is by using the Jacobi equation, which describes 
the displacement £^ between two geodesies under the assumption that £^ 
and ^ are infinitesimal |LC27j . 



Let us show how the Jacobi equation follows from equation (|13p . Using 
Taylor's formula, we can expand each connection coefficients as 

(16) T^(X + = T^(X) + -^(X)C + higher order terms. 

Inserting equation (j 16[) into equation (|13|) and assuming that £ M are 
infinitesimal yields 

(17) e +^-{X)CX y X a + 2T» (T {X)X' J Z V = 0. 

That is, in the above we assume that functions £^ are such that all higher 
order terms , £»£ v , i^t , , CTT 7 . ■■■ can be neglected. Let 

us emphasize that due to this assumption on we can no longer treat 
£ in equation (fTT|) as the exact displacement between two geodesies as in 
equation (fl~2j) . 

Equation (|17p is known as the Jacobi equation for an affine and torsion-free 
connection V. Using the covariant derivative -8- along X, equation (|17j) can 
equivalently be rewritten as 

(18) ^£ + R(£,X)(X) = 0, 
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where R(£,X) the curvature endomorphism of V and £(s) = £ M (s) gfjr|x(s)- 
See for example |Per08| . Prom equations (p~T|) — fjTH]) we see that the assump- 
tions on ^ and £^ have simplified the exact geodesic deviation equation in 
three significant ways: First, unlike the exact geodesic deviation equation, 
equations (fT7|) and ([18]) are linear equations in functions Second, as we 
will see below, equations (|17p and (|18p will be covariant when functions £^ 
transform as components of a tensor. This is a simplification when com- 
pared with the nonlinear transformation rule (|15p for solutions to the exact 
geodesic deviation equation. Lastly, equations (fT7|) and (fT8|) are equations 
along X that only involve evaluations at X(s) on M. 

In the next sections we will study the coordinate invariance of the generalized 
Jacobi equation. As a model for this analysis and to fix notation, let us 
consider in some detail the coordinate invariance for the Jacobi equation. 
Suppose (U,x^) are local coordinates for M and £ = £ M -^jt\ X ( s ) ^ s a curve 
£:/—)• TU along a geodesic X : / — > M. For /i G {0, . . . , n — 1} we define 



(19) « A )^] = e > + ^r^ CT + 2rt t (T ^r. 



That is, the right hand side is the differential operators that appears in the 
Jacobi equation. To simplify the notation we will also write </£/[(£ A )™=o] = 
J}j[£, x ]- We say that a curve £: I — > TU is a Jacobi field inU if tto£: I — > U 
is a geodesic and </^[£ A ] = for all \x G {0, . . . , n — 1}. 



If (U,x^) are overlapping coordinates then 

3x A , 



(20) j£ 



dx^ 



Thus, if we assume that are components for a vector field along X, then 
the definition of a Jacobi field does not depend on the choice of coordinates. 

There are also other ways to derive the Jacobi equation. One approach is 
to start with a geodesic variation around a geodesic X : I — > M. That is, a 
map A: (— ffcb^o) x J — ^ M such that 

(i) s I—?- A(e, s) is a geodesic for each e G (—Eo,Eo), 

(ii) A(0,s) = X(s). 

Then one can show that the tangent of A in the e-direction defines a vector 
field £ : / — > TM, £(s) = d £ A(e, s) | e =o along A, and moreover, the compo- 
nents of £ solves the Jacobi equation (fT7|) . Conversely, any solutions to the 
Jacobi equation on a compact interval, can be written as £(s) = d e A(e, s)| e =o 
for a geodesic variation A. See for example [BD10J. An advantage of this 
derivation is that £ a tangent vector by definition, and the derivation does 
not involve any assumptions like £^ and £^ beeing small or infinitesimal. 



3.3. The generalized Jacobi equation. In the previous section we started 
with the exact geodesic deviation equation (|13|) . did a Taylor expansion of 
T„ a (X + £) (equation (TTo| ) and assumed that £^ and £ M are infinitesimal, 
so that higher order terms f^f, £*f , f^f", ^^? CT , £ M ^ CT , ^"f 0- , • • • can 
be neglected. This gives rise to the Jacobi equation (|17p . The generalized 
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Jacobi equation is derived in the same way, but assuming that only is 
infinitesimal. Under this approximation, equation (|13|) simplifies into the 
generalized Jacobi equation in equation (pQ). 

Suppose (U, x M ) are local coordinates for M, and £°, . . . , £ n : I — > K are 
functions along a geodesic X: I —> U. As for the Jacobi equation we define 
differential operators 

G ^ x ] = e + {^ v x° + rr) + + r) (> + r 

for [i G {0, . . . , n — 1}. We will say that functions £° . . . , ^ n_1 : 7 — y R define 
a generalized Jacobi field in chart U if G^[^ A ] = for \i G {0, . . . , n — 1}. 

Suppose X is a geodesic X : I —> M. Then the Jacobi equation along X is 
locally a linear second order differential equation for components Thus, 
for initial values £ M (so) and £ M (so) there exists a unique Jacobi field £ : 7 — > 
TM along X with these initial values. See for example |Lee97| . In contrast, 
the generalized Jacobi equation Gj/[£ A ] = is locally a non-linear second 
order differential equation for functions The Picard-Lindelof theorem 
implies that for any sq £ I, there exists an neighborhood around sq such 
that the generalized Jacobi is uniquely solvable on this neighborhood from 
initial values £^( s o) and £^(so)- One can bound the size of this neighborhood 
[Har64, Theorem 1.1, Chapter 2]. However, in general there is no guarantee 
that a unique solution exists on the entire interval I as for Jacobi fields. 

For future reference, let us note that differential operators and [£ A ] 

are related by 

(21) Gjj [£ A ] = J£[£ A ] + A^ A ], 



where 

(22) a£k a ] = r(u i r+2^r* i r + ^rrr- 



4. Coordinate invariance of the generalized Jacobi equation 

As described in the introduction, this section contains the main results of 
this note. As in Section [3] we assume that M is a manifold with an affine 
and torsion- free connection V. 

Proposition 4.1. Suppose X : I — > M is a geodesic contained in a chart 
(UyX^) and functions £°, . . . ,C n 1 : I ~ * K define a generalized Jacobi field 
in coordinates x^ . If are coordinates defined by the affine coordinate 
transformation ([7D, then functions ^p;£, X : I ~ * ^ define a generalized Jacobi 
field in coordinates x^ . 

Proof. Let ^ = fpf 7 . Equations ([25]) -([25 jl in Appendix [A] show that 
fpA^ A ] = Ag[£ A ]. By equations flSDJ and we then have fpG^[^ A ] = 
G~[£ A ] and the claim follows. □ 
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Proposition 4.2. Suppose M is a manifold of dimension > 3, X: I — )■ M 
is a geodesic and (U,x^) are local coordinates around X(sq) for some sq G 
/. By shrinking I to a neighborhood of so we can find functions £ M : / — > 
R and overlapping coordinates (U,x^) around X(sq) such that £^ define a 
generalized Jacobi field in chart U , but ffirC does not define a generalized 
Jacobi field in chart U . 



Proof. Let p = X(sq). By proposition 14.11 we may assume that G R n cor- 
responds to p in coordinates x^ . Let (U,x^) be the coordinates determined 
by 



\ 5 - - - 5 y ^ g r P°" ' 



where 



and 5^ and 5^ u are the Kronecker delta symbols. 

Let u = X(sq). Since dimM > 3 we can find vectors v,w G T p M such 
that {u, v, w} is an orthonormal basis with respect to the Euclidean metric 
g = b^ydx^ ® dx v on U. By the Picard-Lindelof theorem we can shrink / 
to a neighborhood of so and find functions £^ : / — >■ R such that define a 
generalized Jacobi field in chart U and 

when it = n^g|r|p, v = v M gfjr| P and u> = w^-£^\ p . Let ^ = $f££ CT . 

Contracting G^[£ A ] = [£ A ] + A^[£ A ] by §§£, applying equations fl2D}-fl2H 
and equations ([23|) ~ (j25|) in Appendix [XI yields 



dx v u*> J m 



dx u dx T dxPdx a 

where all expressions are evaluated at sq. Since Tr pa is symmetric in rpa, 

we have dx ?Q*pdx<* = t t P <t- Moreover, since §§£(p) = 6% and since u,v,w 
are orthonormal it follows that 

G^[£ A ] = -g(2u + w, w) - g(v, w) (2^ + w**) - g(y, 2u + w) wP 



define a generalized Jacobi field in chart U, and the claim follows. □ 



at sq. We have shown that 7^ at sq. Thus functions £^ do not 
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Appendix A. Coordinate transformations 



Suppose X is a curve X : I — > M, and x M and x M are coordinates around 
X(sq) with 9 ®1 X q xP (X(s )) = for some s G If C are functions : J ->■ 

R and functions ^ are denned by £ M = ff^C^, then at X(so) we have 
transformation rules 

(23) e = 

(24) _ pM dx a dx x 

dx^ dvz p _ ar^a^_9x^9P i d 3 ^ 

\ / G~*t d™A 



<9x CT cta; A Sa? 5 <9x A dx p dx u dxPdx x 

These equations show that between coordinates x p and objects 

and -§^r transform as tensors at X{sq) (up to an extra term in the last 
equation). For a general coordinate transformation, these transformation 
rules are considerably more involved. See equation ([2]). 
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